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Abstract 

We study the Euclidean effective action and the full fermion propa- 
gator for a Dirac field in the presence of a scalar field with a domain 
wall defect, in 2 + 1 dimensions. We include quantum effects due to 
both fermion and scalar field fluctuations, in a one-loop approxima- 
tion. The results are interpreted in terms of the quantum stability 
of the zero mode solution. We also study, for this system, the in- 
duced 'inertiaF electric field for the fermions on the defect, due to the 
quantum fluctuations of the scalar field. 
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1 Introduction 



The subject of fermionic fields in the background of defects presents itself in 
many different areas of physics, from textures in superffuid phases of He^ Jl| 
to cosmic strings . Those defects are generally defined as some field configu- 
rations having a nontrivial topological content, and one of the most important 
features of these systems is the existence of localized fermionic zero modes. 
The most interesting property of the zero modes, when transport properties 
are considered, is the fact that they are perfect conductors (namely, they 
have linear dispersion relations for small values of the momenta). 

The most widely known example of this phenomenon is, perhaps, the so 
called Callan-Harvey mechanism || by which a Fermi field (either 'funda- 
mental' or arising as an effective variable in an approximate description of 
a physical system) living in an odd- dimensional spacetime, is coupled to a 
domain wall like defect. This defect is defined as a mass term that changes 
sign on a one-dimensional curve in the constant-time plane. A general result 
for this kind of system is that the fermionic field spectrum has, besides the 
usual continuum part, an extra zero-mass bound-state located precisely on 
the defect. This zero mode corresponds then to a one-dimensional fermion 
which, because of its linear dispersion relation, is capable of carrying electric 
currents. 

When the Fermi field is also minimally coupled to an Abelian gauge field, 
the resulting Chern-Simons current, proportional to the sign of the fermionic 
mass, has an apparent 'anomaly' on the defect, due to the vanishing of the 
mass. This 'bulk' anomaly is, however, exactly compensated by an equal and 
opposite (chiral) anomaly for the fermionic current at the defect, which arises 
because the zero modes are described by a 1 + 1 dimensional chiral theory. 
A similar phenomenon is found at the edge of a droplet of a non-relativistic 
electron gas in two spatial dimensions, in the presence of a magnetic field, 
both in the regime of the fractional and integer quantum Hall effects 0]. Here, 
the Chern-Simons (Hall) current is of course also present, and its apparent 
'bulk' anomaly on the borders is canceled by the gauge anomaly of the chiral 
edge states || ||. Indeed, this anomaly- matching mechanism is usually a 
guiding principle for the construction of the models describing the matter 
fields on the borders, since in two spacetime dimensions the form of the 
anomaly is very closely constrained, though no completely determined, by 
the matter content. 

While for many different physical situations static defects are sufficient 
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to describe the relevant physics, in many others (particularly in condensed 
matter physics) it is of current interest to understand the possible effects due 
to the dynamics of the defects. In particular, it is interesting to investigate 
the interaction of the fermion zero modes with the fluctuating geometry 
induced by the dynamics of the defect. In we generalized the results 
of , which dealt with static defects in 2 + 1 dimensions, to the case of an 
arbitrary, space and time dependent defect, showing that an analog of the 
Callan and Harvey mechanism || still holds. Namely, there exists a chiral 
zero mode, and it is localized, this time on the defect worldsheet. Besides, 
an important new phenomenon is that (even in the absence of any external 
gauge field) if the defect is accelerated, there is an induced fermionic current 
on the defect. This effect is caused by the gravitational chiral anomaly of 
the fermions localized on the defect. This current can be thought of as due 
to the existence of an effective 'inertial' electric field, which can be entirely 
defined in terms of the geometry of the defect, and may also be interpreted 
in terms of the gravitational anomaly for the fermions living on the defect 
worldsheet ||. 

In this paper, we shall consider a related situation: the defect will be 
allowed to fluctuate, but the dynamics of these fluctuations will be, in the 
present case, defined by a scalar field action. To make the system tractable, 
however, we shall deal with configurations where the 'topological charge' of 
the configurations is fixed, corresponding to just one defect. The excursions 
around these configurations will be first taken into account by performing a 
loop expansion around a static defect background. Although the expansion 
itself is a standard procedure, it will demand here the use of a suitable set of 
coordinates, compatible with the (non translation invariant) configuration of 
the system. After showing that the domain wall plus zero mode configuration 
is a consistent classical configuration, we shall show that they are stable under 
quantum fluctuations of the fermionic and scalar fields. 

When the fermion field is a Majorana rather than a Dirac spinor, the 
model becomes supersymmetric. We show that a similar mechanism occurs 
in such a case, although the model itself is less interesting from the phe- 
nomenological point of view, since there is no 'electric' current associated to 
the Majorana fermion. 

The organization of this paper is as follows: in section |2|, we define the 
system and explore some of its fundamental properties. In section || we 
discuss some non-trivial aspects of the loop expansion for the effective action, 
applying this to the derivation of the full fermion propagator including one- 
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loop corrections in section This result is applied to the study of the 
stability of the zero mode solution under the fluctuations. 

The effective inertial electric field, due to the accelerating defect is dis- 
cussed in section j|. 

The supersymmetric version of the model is presented, for the sake of 
completeness, in Appendix A. Section [| presents our conclusions. 

2 The system 

The system we shall deal with may be conveniently defined in terms of an 
Euclidean generating functional: 

ZUiViV) = J T>(pVipVi/j exp[—S((f),ip,i(j) + J d 3 x(j(f> + fjip + iprj)] (1) 

where the (Euclidean) action S(4>, ip, ip) can be written as: 

S{<f>,$,i!>) = S B (<f>) + S F {M,<f>) (2) 

with 

S B (<f>) = Jd 3 x[±(d-^ + V{<P)} (3) 

and 

S F $, = Jd 3 xij{0 + g<j){x)] $ . (4) 

In Equation (||), V(<p) denotes a local potential, which should of course be 
consistent with the existence of a domain wall configuration for 0, and whose 
form will be made more explicit later on. The coupling constant g controls 
the strength of the interaction between the fermions and the scalar field. 

The Euclidean Dirac matrices 7^ have been chosen according to the con- 
vention 

7o = 03 7i = 01 72 = o"2 , (5) 

where the a^s denote the usual Pauli matrices. 

In the absence of external sources, the classical equations of motion for 
this system are: 

[0 + <70(x)ty = (6) 

and 

- d 2 (p(x) + V'(cj)(x)) + gip(x)4>(x) = , (7) 
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where V = ^v. In order to have a fermionic zero mode, one should demand 
the existence of a domain wall like configuration for 0. This, in turn, sug- 
gests that the potential V must be non quadratic, with the simplest choice 
being a quartic function. Thus the domain wall would correspond to an (a?i) 
translation invariant kink solution. This, of course, amounts to an infinite 
energy configuration in the infinite volume case. We have in mind, however, a 
finite system which can nevertheless be regarded as infinite (for calculational 
simplicity) in most of the objects we shall be interested in. 

Let us, for the sake of simplicity, consider a static configuration of this 
type; i.e., (j) = (p(x 2 ). Then (|6]) and (0) reduce to: 

[7o<9 + 7i<9i + 72^2 + 9<P(x 2 M = ( 8 ) 

and 

- + V'tyix*)) + #(z 2 )V(z 2 ) = , (9) 

respectively. To see that a zero mode is indeed a possible consistent solution 
for this coupled system of equations, let us consider a field configuration 
(fi k (x 2 ) verifying 

d 2 if k (x 2 ) 



dx 
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+ V'{ip k {x 2 )) = , (10) 

-2 

with a potential V, such that ipk{x 2 ) is kink-like, changing sign at x 2 = 0, 
say (the 'center' of the kink). Inserting this configuration (p k in the place 
of 4> in @ yields, by the standard Callan-Harvey mechanism, a chiral zero 
mode solution for the fermion field, and this in turn implies the vanishing of 
ifjifj, since this term involves the two chiralities. This means that fk{x 2 ) is 
also an exact solution of the full system of coupled equations (|j) and (|9|). 

Let us now consider the more interesting issue of the quantum dynamics 
corresponding to the above configurations. To that end, we first introduce 
W[j, fj, rj], the generating functional of connected Green's functions, which is 
defined by 

W\j,fj,rj\ = \nZ\j,fj,rj\ (11) 
and the effective action T, the Legendre transform of W: 

r[^,X,x] + W\j,fj,rj[ = J d 3 x[j(x)ip(x) + f)(x)x(x) + x(x)r]{x)] (12) 

where cp, x and \ are the 'classical fields': 

8W , . 5W _ f , 8W 
^ = JJix) X{X) = Wx) t{X) = ~Hx) ( 3) 
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and we adopted the convention that functional derivatives of Grassmann 
objects always act from the left of the corresponding object. 



3 Loop expansion 



A standard one-loop calculation of the effective action T yields [fTO] 

r(y,x,x) = S(ip,x,x) + Fi(<p,x>x) + ••• 



(14) 



where S is the classical action and Ti denotes the order-ft correction, given 
explicitly by 

v x (ip,x,x) = r?(¥>) + r}(v,x,x) (15) 



where: 
and 



-Tr In 



1 1 



6 2 S B (y) 
5(f)(x 1 )S4>(x 2 ) 



-TrlnX? 



with the definition: 



V =@ + g<p 



(16) 



(17) 



(18) 



The two terms contributing to Ti have a different physical interpretation: 
while r™ corresponds to the energy of the 'distorted' Dirac vacuum for the 
fermionic field in the presence of the scalar field ip, T\ takes into account the 
bosonic fluctuations around the classical kink background. We note that, 



besides the standard fluctuation operator 



corresponding to a self- 



&4>(xi)8<i>{x2) 

interacting scalar field, there is also a correction due to the fermion-boson 
interaction. 

As we have already pointed out in the introduction, the field <p will be 
assumed to fluctuate around a background defect. Rather than considering 
arbitrary configurations for that defect, we shall concentrate on those corre- 
sponding to a sufficiently smooth, space dependent step-like defect. More- 
over, the choice of the configuration will also be dictated by the practical 
requirement of making the calculations simpler. Our first simplifying as- 
sumption is to consider a rectilinear defect, with (p being a function of only 
one of the spatial coordinates, say x 2 , and changing sign along the straight 
line x 2 = 0. This will introduce in (|T5|) an explicit dependence on x 2 in 
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both terms. This dependence is obviously carried by the operator T> defined 
in (0) but also an x 2 dependence will be induced in the double functional 
derivative of Sb appearing in T\, since it has to be evaluated on the classical 
configuration for ip. To simplify the task, taking advantage of the translation 
invariance on the x$x\ plane, one should be able to disentangle the dynam- 
ics of the fermions in two pieces: one depending only on x 2 , and the other 
having support exclusively in the plane x 2 = 0. To this end, we note that 
the operator D may be rewritten as: 

V = (a+ 0)V L + (a f + 0)V R (19) 

where a, are operators acting on functions of x 2 , 

a = d 2 + gip(x 2 ) , a f = -d 2 + g<p{x 2 ) , (20) 

Vl,r are projectors along the eigenspaces of the 'chirality' matrix j 2 : 

V L = ^(1 + 72) V R = i(l- 72 ), (21) 

and is the two-dimensional Euclidean Dirac operator corresponding to the 
two coordinates xq and x%, which we denote collectively by x, namely 

I = 1 ■ d = la d a , a = 0, 1 . (22) 

Expression ( |T9"D suggests that one could get rid of the x 2 dependence 
of the fields by a suitable expansion in the modes of some operator, thus 
obtaining a 'dimensional reduction' from the three dimensional spacetime to 
the two dimensional one described by x. As D itself is not Hermitian, we 
may use instead any of the positive Hermitian operators 7i or TC: 

n = v ] v , n = vv ] (23) 

to define the modes. 

Using the expression (|19| ) for T>, leads to the explicit form for TC and 7i: 

n = -d 2 + <?y - g l2 d 2lf , n = -a 2 + <?V + rt&v ■ (24) 

These expressions enable us to write the contribution to the effective 
action either in terms of TC or TC: 

r$((p) = -^TilnTC = -^TrlnW. (25) 
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Taking the trace over Dirac indices, and using any of the two previous rep- 
resentations, we end up with: 

r?(¥>) = -\^H-d 2 + g 2 V 2 + gd 2V )- l -Ti\n{-d 2 + g 2 V 2 -gd 2V ) (26) 

where the trace is now meant only on functional space. 

On the other hand, in T\ we may factorize a bosonic fluctuation determi- 
nant, so that 



r? = -Trln 
1 2 



5 2 S B 



6(p(x 1 )Sip(x 2 ) 



+ -Trln 

2 



5(x 1 -x 2 )-2g A ip (x 1 ,y)x(y)V (y, x 2 )x(x 2 ) 



(27) 



where: 



S 2 S 
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6<p(y)8(p(x 2 ) Jy 5(p(x 1 )5ip(y) 



5 2 S f 



-A v (y,x 2 ) = 8(xi -x 2 ) . (28) 



With equation fl2"7p in mind, we note that the one-loop effective action, 
r 1; may also be written as: 



r (o) + r (D ^ 



(29) 



where represents the contribution that would correspond to a system 
consisting of a free fermion in a non-trivial ip background plus the vacuum 
energy due to the self-interacting scalar field (p, with no term taking into 
account the interaction energy between scalar and fermion fields: 



(0) 



- - Tr In 

2 



— Tr In 

2 



+ - Tr In 

2 



d 2 -d 2 + g 2 (p 2 (x 2 ) - gd 2 ip(x 2 ] 

d 2 -d 2 + g 2 ip 2 (x 2 ) + gd 2 <p(x 2/ 
~d 2 - d\ + V"{ V {x 2 ))~ 



(30) 



The term, instead, is a measure of the interaction energy between bosons 
and fermions, vanishing when g — ► 0: 



-Tr In 

2 



8(x 1 -x 2 )-2g A tp (x 1 ,y)x(y)T> (y,x 2 )x(x 2 ) 



■ (31) 
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In what follows, we shall discuss the contributions and T^' sepa- 
rately. We first make the observation that the 'vacuum' term may be 
conveniently put in the following form: 



(0) 



— Tr In 

2 



-d 2 + h f 



Tr In 

2 



-d 2 + h f 



+ ^Trln 

2 



-d 2 + h v , (32) 



where the dependence on the domain wall profile (fi(x 2 ) and on the potential 
V(ip(x 2 )) is encoded in the operators hf, hf and h v , which are defined by 



hf = a)a , hf = aa) , h v = —d 2 + V"(ip(x 2 )) 



(33) 



These operators play the role of one dimensional quantum mechanical 'Hamil- 
tonians', whose eigenvalues appear as parameters in the corresponding func- 
tional determinants. The eigenvalues and eigenvectors of these operators are 
defined by the equations 

hf^ n {x 2 ) = \ 2 n i) n (x 2 ) h f i) n (x 2 ) = \ 2 n 4) n (x 2 ) 

K^ n {x 2 ) = H 2 n £ n {x 2 ) , (34) 

where all the eigenvalues are written as squares (of real numbers) to empha- 
size the fact that the corresponding operators are positive. 

We now start to consider some simplifying assumptions regarding the 
shape of the defect. Noting that 



r?(y?) = ~TY\n(-d 2 + h) -\Tr\n(-d 2 + ~h) 

Z Z 



(35) 



where h = a' a, h = aa^ and o = d 2 + gtp, we assume that only h has a 
zero mode |0), which of course verifies the first order equation a\0) = 0. The 
rest of the h and h spectra coincide, and we demand it to be a continuum, 
separated by a finite gap k 2 > from the zero mode. The resulting equations 
uniquely fix the form of ip to be: 



K 



9 



tanh(/-€£ 2 ) , 



(36) 



where the constant k remains arbitrary. The first order equation for the zero 
mode can be integrated by quadratures, yielding: 



^0(^2) 



£2 



^ [cosh(ra 2 )] 1 , 



(37) 
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which has been normalized to (0 1 0) = 1. For the previous ip field profile, we 
see that the explicit forms of the operators hf and hf become 

h f = -d\ + K 2 [2tanh 2 (ra 2 ) - 1] h f = -d\ + k 2 . (38) 

Although it is not immediately evident from (|38"D, the spectra of hf and 
hf coincide, except for the zero in hf. Indeed, any non-zero mode of hf 
may be used to generate an eigenvector of hf (with identical eigenvalue), 
by application of a\ and reciprocally, applying a to a non-zero mode of hf 
produces an eigenvector of hf. On the other hand, by assumption, only hf 
has a zero mode. Taking this into account, we see that may be written as 

r« = -Tr 2+1 ln(-a 2 - dl + k 2 ) - iTr 1+1 \n(-d 2 ) (39) 

where the dimension of the spacetime where the functional trace is defined 
has been explicitly shown. Of course, we may also go back to a 'spinorial' 
representation: 

r; = -Tr 2+1 ln(^+ K ) - Tr 1+1 \n(fiV L ) , (40) 

where the zero mode contribution has been singled out, and the trace affects 
spinorial indices also. 

Fixing (|36|) to be the form of tp, we now need to make sure that it is a 
minimum of Sb- This may be achieved by a judicious choice of the potential 
V. As already advanced, the simplest choice is to use a quartic potential; 
indeed, assuming the form of V(ip) to be: 

v{r) = \i* 2 - , (4i) 

then there is a classical kink solution 

(p(x2) = jtt tanh(vA (1x2) ■ (42) 

In order to match the configuration (|36| ) (something we require in order to 
simplify the spectra of h and h), we should demand: 

fi = — , VXfi = k (43) 
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so that the scalar potential is: 



V{<P) = y (V " (-/) ■ (44) 



With V as given by (0), we see that h 9 is 

h v = -d\ + 6K 2 tanh 2 (ra 2 ) -2k 2 . (45) 

The spectrum of h v can also be found by algebraic means. We note that h v 
can be factorized similarly to hf. 

h v = tfb (46) 

where we introduced: 

b = 82 + 2Ktanh(ra2) , b = — 82 + 2Ktanh(ftX2) • (47) 

This already shows that there will be a zero mode £0 for h v , which verifies 
^o = 0, and is given by: 

I~3k 

CoW = Jj[cosh(ra 2 )]- 2 , (48) 

when normalized to (£o|£o) = 1- Regarding the spectrum above the zero 
mode, we may take advantage of the fact that the operator h v = btf verifies: 

h v = h f + 3k 2 (49) 

to infer that h v will also have another bound state (with eigenvalue 3k 2 ), 
and then a continuum starting at 4k 2 . The form of these extra states may of 
course be found by applying the operator to the corresponding eigenstates 
of hf. 

Thus the spectrum of h v contains two discrete states of eigenvalues and 
3k 2 , plus a continuum starting at 4k 2 . Using this information in combination 
with the already known form of T", we see that 

ri 0) = -iTr 1+1 ln(-«9 2 ) + iTr 1+1 ln(-«9 2 ) + ~Tr 1+1 ln(-d 2 + 3k 2 ) 

- Tr 2+1 ln(-d 2 - d\ + k 2 ) + ^Tr 2+1 ln(-d 2 - d\ + 4k 2 ) , (50) 

where all the terms traces affect only functional space. Note that the first 
two terms, which correspond to the fermionic and bosonic zero modes, cancel 
out exactly. 



11 



4 The full fermion propagator 



The full fermion propagator Sp may be determined by evaluating the inverse 
of the second functional derivative of the (one-loop corrected) effective action, 
evaluated at zero fermionic field, namely, 



S F = V- 1 (51) 

where 

T)(t.ii\ = - 



V^iV) = ? J x=x=o = V(x,y) + V 1 (x,y) (52) 



with 

V 1 (x,y) = g 2 A^y,x)V-\x,y). (53) 

It is evident that V plays the role of a 'corrected Dirac operator', since it in- 
cludes a one-loop contribution. Contrary to what happens for the translation 
invariant case, we see that the use of a standard Fourier representation for all 
the coordinates does not simplify the treatment. Rather, a mixed expansion 
using plane waves for the x coordinates and suitable eigenfunctions for the 
dependence on x 2 is more adequate. Indeed, in terms of these eigenfunctions 
and eigenvalues, we may formulate the expansion for the propagator V^ 1 as 
follows 

+ ^[(x 2 \ij n )(i) n \y 2 )-^-V R + {x 2 \^ n )^ n \y 2 )^-^-V L 



n=l 

00 A A 1 

+ EK^n) (^2)^72^+ {x 2 \^ n )^ n \y 2 )-^-V L ] . (54) 

n=l P "r" *n V + *n ) 

We are using a discrete sum notation here for the sake of simplicity; it should 
be noted, however, that with our assumptions on the potential, there is only 
one discrete state, ip . We shall take care of this, using a more explicit 
notation when necessary. 

On the other hand, for we use instead the expansion: 

r c] 2 % 00 1 

kA^\v) = -j (2^)2 e*-<*-» E<^IQ<^>^^ ■ (55) 
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Of course, the same remark on the number of discrete modes (now there are 
two of them), applies here. 

In terms of the kernels defined above, we see that 

V ^ x ^ = - jrf^h s = 9 2 \{y,*)(x\v- l \y) , (56) 

whose explicit form in the mixed Fourier representation becomes: 

Vx{x,y) = J-^e ik <*-*h k (x 2 ,y 2 ) (57) 

with 

{oo 
Y ^(^2)^0(^2)^(2/2)^0(2/2) J {k; fin, 0)V R 
n=0 

00 

+ Y [^(^2)^(^2)^(2/2)^(2/2) J(k fin, ^m)V R ] 
n=0,m=l 

00 

+ Y €nM$ m fa)t* (2/2)^(2/2) j(k\fl n ,\ m )V L 
n=0,m=l 

00 

n=0,m=l 

00 1 

+ Y ^(2:2)^(^2)^(2/2)^(2/2) I{k; fi n Am)V L > , (58) 

n=0,m=l \ 

where we used a discrete notation for the sum over eigenvalues, but of course 
an integral over the continuous part of the spectrum is implicitly assumed. 
More precisely, we only have one discrete zero mode (Ao = 0) for the fermionic 
part, plus two discrete modes with //jj = and fi\ = 3k 2 for the bosonic part 
of the sums. 

The functions I and J appearing in ([58]) may be interpreted as 1 + 1 
dimensional loop integrals. They are explicitly given by 

up m2 

(2^) 2 [(k-p) 2 + M?](p* + M$) 
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and 



J{k;M 1 ,M 2 



d 2 p —i § 

(2vr) 2 [(k - p) 2 + M 2 ](p 2 + M 2 



(60) 



These (UV convergent) functions may be evaluated by using standard meth- 
ods, what yields 



I{k;M 1 ,M 2 ) 



In 



Airy (k 2 + AM 2 )(k 2 + A 2 ) 



(k 2 + AM 2 )(k 2 + A 2 ) 



and 

J{k;M 1 ,M 2 ) 
where 

a' 



k 2 + M 2 + M 2 + \f(k 2 + 4M 2 )(F + A 2 ) 
k 2 + Ml + Mf - 

i % 



Ank 2 (a~ 



a + ln( 



a + — 1 , « — 1 . 

I — a m(- 



af 



of 



1 

2k 2 



k + 2MA =f \/(k 2 + 4M 2 )(k 2 + A 2 ) 



(61) 



(62) 



(63) 



and we used the notation: M 



M1+M2 



2 , A = Mi — M 2 . 
We now consider the effect of including this correction into the Dirac 
operator. The corrected Dirac operator D, equation (^), in a representation 
where the first two components have been Fourier transformed, will have a 
kernel 

(x 2 \V\y 2 ) = (x 2 \(a + ij£)\y 2 )V L + (x 2 \(J + ify\y 2 )V R + -y K (x 2 ,y 2 ) (64) 

where the 'free' part is of course local and translation invariant, since 

(x 2 \a\y 2 ) = a(x 2 )5(x 2 - y 2 ) (65) 

and the analogous relation for a'. The correction 7^, in turn, will be non 
local, with an structure which may be described as follows: 

lk{x 2 ,y 2 ) = u k (x 2l y 2 )V L + u k (x 2 ,y 2 )V R 

+ v k (x 2 ,y 2 )iftV L +v k (x 2 ,y 2 )iftV R , (66) 

where the functions Uk,Uk,Vk,Vk have a rather complicated expression, as 
may be inferred from (j58|). 
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To study the stability of the chiral zero mode solution, we note that the 
zero mode |^o) °f the 'free' Dirac operator, £>, may be written as: j^o) — 
VLipo(x2)\k) , were \k) is a zero mass spinor, solution of the 1 + 1 dimensional 
Dirac equation. If we apply now the corrected Dirac operator D to the same 
state l^o), we see that: 

r+oo 

{x 2 \V\* ) = 0+ / dy 2lk (x 2 ,y 2 )(y 2 \Vo) , (67) 

and let us now quantify the effect of the rhs. From the explicit form of 7^, 
we see that only the terms proportional to Vl may contribute. Moreover, 
the integrals over the y 2 coordinate will be of the form: 

/+00 „ 
^2/2 ^(2/2)^(2/2)^0(2/2) = turn- (68) 
-co 

The explicit forms of the functions appearing in the integral are 

^0(2/2) = ^sech(/«/ 2 ) (69) 

for ip , and 

#U(wO = -^jV^y* (70) 



2tt 

since hf is a (shifted) free particle one dimensional Hamiltonian. Of course, 
> k 2 . Regarding £„, it has different profiles depending on whether n — 0, 
n = 1 or belongs to the continuum: 

I~3k 

6(2/2) = W — sech 2 (/q/ 2 ) , (71) 
£1(2/2) = bH {y 2 ) = iV(-a 2 + 2Ktanh(«y 2 ))cosh(^ 2 ) = ^ 

1?2) 

and for states in the continumm, 

£t(y 2 ) = -L (-i^2-4K2 + 2«tanh(« lfe )) e -*VS=^i« , (73) 

V 



where /u 2 > 4k 



2 
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We now show that £„ jm is exponentially small. When n = 0, we have: 

dy 2 £0(2/2)^(2/2)^0(2/2) = 



f \/f l ^ C0S [V^ ~~ « 2 ?/ 2 ]sech 3 (/q/ 2 



3\/7r ! X m , 2 




(-)^ech[-W(^-l] (74) 



which is exponentially decreasing with A m . For the n = 1 case, we have 
instead: 

dy2 £1(2/2)^(2/2)^0(2/2; 




2 - 1)sech[ iv ( v )2 - 1] (75) 

which is also decreases exponentially. 

When n belongs to the continuum, a very similar inequality shows that 
£„ im is also exponentially small, and indeed much smaller than in the n = 
case. 

All this shows that the free zero mode l^o) is still a zero mode of the cor- 
rected Dirac operator, modulo exponentially small corrections. It is notewor- 
thy that the parameter that regulates the rate of decrease of the exponentials 
is k. Thus, the sharper the domain wall, the more negligible the departures 
shall be. Of course, the fact that the correction is so small suggests the use 
of the free zero mode as a starting point for a perturbative evaluation of the 
corrected one, by repeated application of the corrected Dirac operator. 

It is interesting to see what happens when the gap between the zero modes 
and the higher modes is very large, at the level of the interaction part of the 
effective action I 1 ] 1 , equation fl3ll). When k becomes very large, we may 
approximate the propagators as follows: 

V(x,y) ~ ^0(^2)^0(2/2)^ (x,y) , V(x,y) = -(x\p V R \y) 

\(x,y) ~ £0(^2) £0(2/2) \(x, y) , A v (x,y) = - (x\d~ 2 \y) . (76) 

Also, by a similar argument to the one used for the fermionic propagator, 
we see that the 'classical fields' x appearing in (^l|) may be replaced by the 
projections: 

X(x) ~ ^ (x 2 )x(x) (77) 
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where x(x) is an arbitrary two dimensional spinor. Using these approxima- 
tions, we find that 



lf } Trln 

1 2 



8(x -y)~ 2g 2 ^j -£ (x 2 )<j)l(x2) G(x, y) 



where 



On the other hand, from the relation between i/jq and £o, we note that 



G(x,y) = / A v (x,z)x(z)V(z,y)x{ 



^0(^2)^0(^2) 
and finally, this implies: 



g(^2|Co)(Co|y2> , 



rS 7) Trln 

1 2 



5(x-y)-2g\-G(x,y) 



(78) 
(79) 

(80) 
(81) 



where the x 2 dependence has disappeared. 



5 Induced electric field 

Let us discuss here the application of some general results [0] regarding the 
induced electric field seen by the fermionic zero mode, to show that, for this 
kind of configuration, and in this approximation, that electric field vanishes 
on the average. 

The effective electric field will only be non-zero when the defect is ac- 
celerated, and this can only happen, in our situation, when we include the 
fluctuations in the scalar field. We recall that, when the fluctuations are small 
and the defect is rectilinear, the formula that yields the effective electric field 
E e ff for the fermionic zero mode in a given scalar field configuration is: 



1 r 



d^dfr] - (dodtrjY (82) 



where 7](xq,Xi) is the fluctuation of tp around its static kink configuration 
with center in x 2 = 0. Expression (|82|) is the value of the electric field for 
a given tp configuration; it is of course interesting to consider its quantum 
average (E e ff), the average defined by the Euclidean action of the model. 
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Before considering the specific form of the average, we note that the integral 
of E e ff necessarily vanishes, since 



dxqdxi E, 



eff 



dx dxi d^rjO^r] — (d dirjy 







(83) 



by an integration by parts. This has consequences for the spacetime av- 
erage of (E e ff), but it leaves room for the existence of a locally non-zero 
(E e ff). The form of this local effective electric field can be written, in our 
approximation, as follows: 



(£ e// (x ,zi)) 



8 2 d 2 



d 2 



d 2 



8xq dy\ dx dx\ dyodyi 



\{xo,xi;yo,yi)\ 



(84) 

Now it is evident, 
4|) vanishes. Indeed, 



where A v is the 'fluctuation operator' defined in 
since is translation invariant in xq and x\, that 
translation invariance implies that A v can only depend on the differences of 
its arguments, and this makes it possible to show that the two terms on the 
right hand side of flSif ) cancel out. 

This negative answer for the quantum average might of course be expected 
on symmetry grounds. However, it can be used the other way around, to look 
for the existence of a non-vanishing average by considering different scalar 
field configurations violating the invariance under translations in more than 
one coordinate. 



6 Conclusions 

We have studied a system consisting of a Dirac fermion interacting with a 
real scalar field, to understand the effect of quantum fluctuations on the zero 
modes due to the Callan and Harvey mechanism, to one-loop order. We have 
shown that the effect of the fluctuations on the zero mode is exponentially 
suppressed, and hence negligible, what signals the perturbative stability of 
the classical zero mode solution under quantum fluctuations. 

The effective inertial electric field induced on the fermions vanishes to 
this order; and this also points to the direction of stability, since an strong 
electric field on the zero mode may have meant that an extra interaction 
should have been included. 

Based on these results, we conclude that instabilities may very likely arise 
from different configurations, possible with more than one defect, and with 
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translation invariance along only one spacetime coordinate. Also, it remains 
an open question whether a nonperturbative treatment for the quauntum 
coupled system is possible. Nonperturbative results do indeed exist for this 
system, although they deal with the external field problem only . 

Results on these systems, as well as on numerical nonperturbative solu- 
tions of the coupled system, will be reported elsewhere. 
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Appendix A: Supersymmetric model 

We discuss here some particularities of the (simplest) supersymmetric version 
of the model discussed previously. Taking into account the fact that the 
bosonic part of the model has to be a real scalar Q and we want a non- 
trivial interaction between boson and fermion fields, we are lead to consider 
a 2 + 1 dimensional Wess-Zumino model supplemented by a mass term and 
an interaction term. 



Following the notation and conventions of [12], we consider the action 



S — Swz + S m + Sg 
where Swz is the free, massless, Wess-Zumino action 



S 



wz 



d x 



+ -F 2 
2 



(85) 



(86) 



where we used a Minkowski space metric r\^ v = diag(— 1, 1, 1), and ip is a 
(real) Majorana fermion. The actions S m and S g are mass and interaction 
terms, respectively: 



Sn 



d 3 x — ( Fcj) i/ji/j 



This must be so in order to have domain wall like solutions and fermionic zero modes. 
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So 



d 3 x g (F(J) 2 - ijjijj<j) 



(87) 



As usual, the role of the auxiliary field F is to make the action S of-shell 
invariant under the supersymmetry transformations, which, as shown in |12| , 
are given by: 



50 
5F 

dtp 



e & 



be + Fe 



Supersymmetry is preserved, although now on-shell, if the auxiliary field F 
is 'integrated out'. This yields the equivalent action 



S 



I 



d 6 x 



1 K 

2^2 ( 



+ 



/.2\2 



(89) 

In order to compare with the non supersymmetric action, it is convenient 
to make the field redefinition 4>{x) — > <£>(x), with <p{x) = <fi(x) + so that, 
in the new variables, the action is: 



S 



d 3 x 



1 



1 -r 



—d^ip -^W- j(v 2 - (-) 2 ) 2 - 9VW 



(90) 



This action is almost identical to its non-supersymmetric version, differing 
only in their fermionic sectors. 
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